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Abstract 

The changes of dynamical behaviour of a single fullerene molecule inside an armchair carbon nanotube caused 
by the structural Peierls transition in the nanotube are considered. The structures of the smallest C20 and Fe@C2o 
fullerenes are computed using the spin-polarized density functional theory. Significant changes of the barriers 
for motion along the nanotube axis and rotation of these fullerenes inside the (8,8) nanotube are found at the 
Peierls transition. It is shown that the coefficients of translational and rotational diffusions of these fullerenes inside 
the nanotube change by several orders of magnitude. The possibility of inverse orientational melting, i.e. with a 
decrease of temperature, for the systems under consideration is predicted. 



1 Introduction 

The structure and elastic properties of carbon nanotubes are 
studied in connection with the perspectives of their applica- 
tions in nanoelectronic and nanoelectromechanical devices 
and composite materials, and are also of fundamental inter- 
est, particularly for physics of phase transitions. For exam- 
ple, superconductivity (Tj, commensurate-incommensurate 
phase transition in double-walled nanotubes [2 1, spontaneous 
symmetry breaking with formation of corrugations along 
nanotube axis |3| and structural Peierls transition in arm- 
chair nanotubes [4-9 1 have been considered. In the present 
Letter, we consider a fundamentally new phenomenon re- 
lated to phase transitions in nanosystems. In other words, 
we consider the possibility of inverse orientational melting 
for molecules encapsulated inside nanotubes caused by struc- 
tural Peierls transition in the nanotubes. 



The possibility of Peierls transition in carbon nanotubes 
was first considered in [4j. As a result of this transition, arm- 
chair nanotubes become semiconducting at low temperature, 
and Peierls distortions lead to the Kekule structure (see Fig- 
ure [TJ with two essentially different C-C bond lengths and 
a triple translational period (three times more hexagons in 
the translational unit cell). In previous studies, the Peierls 
gap fl5}f7) and the temperature of the transition to the metal- 
lic phase with equal C-C bond lengths [4, 5 8] were esti- 
mated. Recently, the Kekule structure was calculated for 
the ground state of an infinite armchair (5,5) nanotube by 
PM3 semiempirical molecular orbital calculations [9 1. It was 
shown that, for the (5,5) nanotube, the difference between C- 
C bond lengths for semiconducting phase is 0.03 A, whereas 
the difference between nonequivalent C-C bond lengths for 
metallic phase is only 0.006 A |9|. Note that density func- 
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Figure 1 - Calculated Kekule structure corresponding to the 
ground state of the (8,8) nanotube. Nanotube axis is shown by 
the dashed line. 

tional theory (DFT) calculations for the (5,5) nanotube of a 
finite length 110111 II also gave a 60 atom periodicity of phys- 
ical properties on the length of nanotube segment which is 
consistent with Kekule structure for infinite armchair nanotu- 
bes. Moreover, X-ray crystallographic analysis of chemically 
synthesized short (5,5) nanotubes [12J shows the Kekule 
bond length alternation pattern, which was in good agree- 
ment with DFT and PM3 calculations also performed in lfl2l . 
By the example of the infinite (5,5) nanotube, it was demon- 
strated that the structural Peierls transition connected with 
spontaneous symmetry breaking takes place not only with an 
increase of temperature, but also can be controlled by uniax- 
ial deformation of armchair nanotubes J5J. 

A dynamical behaviour of molecules encapsulated inside 
nanotubes can correspond to the following regimes: oscilla- 
tions about a fixed position and/or a fixed orientation of the 
molecule (regime A), hindered motion along the nanotube 
axis and/or rotation of the molecule (regime B) and free mo- 
tion and/or rotation of the molecule (regime C). In the present 
Letter, we show the possibility of changes of the dynamical 
behaviour of molecules encapsulated inside armchair nanotu- 
bes as a result of the Peierls transition in the nanotube struc- 
ture. In other words, these changes can include switching be- 
tween the regimes A and B, switching between the regimes B 
and C, and the changes in diffusion coefficients correspond- 
ing to the hindered motion and/or rotation of the molecule 
(regime B). The considered changes are possible in the case 
where the regime B takes place for at least one phase of the 
nanotube, i.e. the temperature 7p of the Peierls transition 
should correspond to the temperature range of the regime B 
(the hindered motion and/or rotation of the molecule). In 
other words, the temperature Tp should be of the same order 
of magnitude (or a few orders of magnitude less) as energy 
barriers AE for motion and/or rotation of the molecule in- 
side the nanotube at this phase. Note that inverse melting 
of motion and/or rotation of the molecule is possible, if the 
Peierls transition from the high- to low-temperature phase of 



the nanotube occurs with switching from the regime B to the 
regime C or switching from the regime A to the regime B. 

Estimations showed that the Peierls transition temperature 
is T P ~ 1-15 K dSH. According to calculations lfl3lfT4l . 
the barriers of the value close to this temperature range 
were obtained for rotation of the fullerene Cgo inside the 
C6o@C24o nanoparticle. It was also found that the changes 
of bond lengths of the fullerene C240, the outer shell of these 
nanoparticles, within 0.06 A lead to an increase of the barri- 
ers for rotation by more than an order of magnitude (11311141 . 
The changes of the nanotube bond lengths caused by the 
Peierls distortions are of the same order of magnitude (about 
0.02 A for the (5,5) nanotube 0). Note also that the size 
of an encapsulated molecule, and therefore, the nanotube ra- 
dius cannot be too large, since the magnitude of the Peierls 
distortions decreases with an increase of the armchair nano- 
tube radius 1 6 1 . 

Thus, taking into account the above considerations, we 
have chosen the smallest fullerene C20 and the magnetic 
endofullerene Fe@C2o to investigate changes in the dynam- 
ical behaviour of molecules inside nanotubes at the Peierls 
transition. It has been shown that the (8,8) nanotube is the 
smallest armchair carbon nanotube which can encapsulate 
the fullerene C20 132) • A carbon nanotube with the fullerene 
C20 inside was also used as a model system to simulate a 
drug delivery via the nanotube [16|. 

This Letter is organized as follows: Section 2 presents the 
DFT calculations of the structure of the C20 and Fe@C2o 
fullerenes and the PM3 calculations of the structure of the 
(8,8) nanotube. Section 3 presents the semiempirical calcula- 
tions of the barriers for motion and rotation of the fullerenes 
inside the nanotube. Section 4 is devoted to the dynamical 
behaviour of molecules inside the nanotubes. Our conclu- 
sions are summarized in Section 5. 



2 Fullerene and nanotube structures 

Structures of the C20 and Fe@C2o fullerenes have been cal- 
culated using the spin-polarized density functional theory 
implemented in NWChem 4.5 code ifFTl with the Becke- 
Lee- Yang-Parr exchange-correlation functional (B3LYP) 
H18II19I . Eighteen inner electrons of the iron atom are emu- 
lated with the help of the effective core potential — CRENBS 
ECP [20 1 (only 8 valence s-d electrons are taken into ac- 
count explicitly). The 6-3 1G* basis set is used for describing 
electrons of the carbon atoms. 

The ground state of the fullerene C20 is found to be a sin- 
glet state and has £>2/, symmetry. The calculated energy of 
the triplet state of the fullerene C20 is found to be 64 meV 
greater than the energy of the ground state. The ground state 
of the endofullerene Fe@C2o is found to be a septet state and 
has C211 symmetry. The calculated structures of the ground 
states of the C20 and Fe@C2o fullerenes are shown in Fig- 
ure |2] The iron atom locates in the centre of the endofulle- 
rene. The smallest and the greatest distances between the 
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Figure 2 - Calculated structure of the C20 (a) and Fe@C2o (b) fullerenes. The bond lengths are in angstroms. The distances between 
the carbon atoms and the fullerene centre (in angstroms) are denoted in italics. The atoms which have smaller and greater distances to the 
fullerene centre are shown by the open and filled circles, respectively. C2 symmetry axis is shown by the dotted line. 



carbon atoms and the C20 fullerene centre increase by 6 and 
1%, respectively, as a result of the iron atom encapsulation. 

The semiempirical method of molecular orbitals modified 
for one-dimensional periodic structures lETI with PM3 pa- 
rameterization [22] of the Hamiltonian has been used to cal- 
culate the structure of the (8,8) nanotube. The method was 
used previously for calculating the Kekule structure of the 
(5,5) nanotube ground state and for studying structural transi- 
tions controlled by uniaxial deformation of this nanotube ||9). 
The adequacy of the PM3 parameterization of the Hamilto- 
nian has been demonstrated ll23l by the calculation of bond 
lengths of the Cfio fullerene with J/, symmetry: the calculated 
values of the bond lengths agree with the measured ones 11241 
at the level of experimental accuracy of 10~ 3 A. The calcu- 
lated Kekule structure of the (8,8) nanotube ground state is 
shown in Figure Q] The difference between the lengths of 
short and long bonds of this Kekule structure of the (8,8) 
nanotube is close to such a difference of the (5,5) nano- 
tube H). The Peierls distortions include also radial distor- 
tions of the armchair carbon nanotube with periodicity of 
half of the translational period of the nanotube (for details 
see |9]])- In the case of the (8,8) nanotube, the longest nano- 
tube radius is 0.547 nm, while the shortest radius is 0.544 nm. 



3 Fullerene-nanotube interaction 

The structures of the fullerenes C20 and Fe@C2o obtained in 
section 2 have been used for finding the ground state posi- 
tion and for studying motion and rotation of these fullerenes 
inside the (8,8) nanotube. Two structures of the (8,8) nano- 
tube have been considered: the Kekule structure calculated in 



section 2, and the structure with all equal bonds 1 .423 A in 
length (so that to be equal to the average bond length of the 
calculated the Kekule structure). 

The calculations of the interaction energy between walls 
of double-walled carbon nanotubes showed that the barriers 
for motion of the short wall relative to the long wall are very 
sensitive to the length of the long wall 11251 . Thus, the size of 
the system is too large for ab initio calculations. The analo- 
gous problem exists also for the considered case of a single 
fullerene inside a nanotube. Therefore, the interaction be- 
tween carbon atoms of the fullerenes and the nanotube at the 
interatomic distance r is described by the Lennard- Jones 12- 
6 potential 



(1) 



with the parameters e = 2.755 meV, o = 3.452 A. These pa- 
rameters of the Lennard- Jones potential for the fullerene- 
nanotube interaction are obtained as the average values of 
the parameters [26] for fullerene-fullerene and fullerene - 
graphene interactions, in accordance with the procedure de- 
scribed in [26]. Here, the Lennard- Jones potential is used 
for calculating the potential surface of the interaction energy 
£\v between the fullerene and the infinite nanotube, and we 
believe that this gives adequate qualitative characteristics of 
the potential surface shape. The cut-off distance, r = r c of 
the Lennard-Jones potential is taken equal to r c = 15 A. For 
this cut-off distance the errors of calculation of the interac- 
tion energy £\y between the fullerenes and the (8,8) nano- 
tube and the barriers for relative motion and rotation of the 
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Figure 3 - The interaction energy £\y (in meV) as a function of the fullerene C20 displacement z along the nanotube axis and the 
angle (p of the fullerene rotation about the nanotube axis, (a) The (8,8) carbon nanotube with the Kekule structure; (b) the (8,8) carbon 
nanotube with the structure of metallic phase. The energy is given relative to the energy minima. The equipotential lines are drawn at an 
interval 0.2 meV. 



fullerenes inside the nanotube are less than 0.1%. Both the 
fullerenes and the nanotube are considered to be rigid. An ac- 
count of structure deformation is not essential for the shape 
of the potential surface both for the interwall interaction of 
carbon nano tubes 11251 1271 and the intershell interaction of 
carbon nanoparticles [13. 14]. For example, the account of 
the structure deformation of the shells of C60 @C24o nanopar- 
ticle gives rise to changes of the barriers for relative rotation 
of the shells which are less than 1% 113111411 . It should also 
be noted that the symmetry of interaction energy as a func- 
tion of coordinates describing relative positions of interact- 
ing objects is determined unambiguously by symmetries of 
the isolated objects and does not change if the symmetries of 
the objects are broken because of their interactions. 

The ground state interaction energies between the C20 
and Fe@C2o fullerenes, and the (8,8) nanotube with Kekule 
structure are found to be —1.596 and —1.598 eV, respec- 
tively. The angles between the C2 symmetry axes of the C20 
and Fe@C2o fullerenes and the nanotube axis at the ground 
states are49.6° and 53.1°, respectively. The metastable states 
with the C2 symmetry axes of the fullerenes perpendicular to 
the nanotube axis are also found for both C20 and Fe@C2o- 
At the metastable states, the interaction energies are greater 
by 4.58 and 3.04 meV than the ground state energies, for C20 
and Fe@C2o, respectively. 

The potential surfaces of the interaction energy between 
the C20 fullerene and the nanotube, Ey?((p,z), as functions of 
the relative displacement of the fullerene along the axis of the 
nanotube z and the angle of relative rotation of the fullerene 
about the axis of the nanotube (p are presented in Figure [3] 
for both the considered structures of the (8,8) nanotube. In 
the general case, diffusion of a fullerene along the nanotube 
axis is accompanied by rotation of the fullerene. Our calcu- 
lations show that the barriers for rotation of both fullerenes 



about the axes which are perpendicular to the nanotube axis 
lie between 3 and 23 meV for any orientation of the fullerene 
and for both the considered structures of the (8,8) nano- 
tube. These barriers are significantly greater than the bar- 
riers for rotation of the fullerenes about the axis of the nano- 
tube (shown in Figure [3] for the C20 fullerene). Therefore, 
diffusion of the fullerenes along the nanotube axis is accom- 
panied only by rotation about the nanotube axis. Thus, the 
minimal barrier AE^ for diffusion of the fullerenes along 
the nanotube axis is the barrier between adjacent minima 
of the potential surface £\y(<P,z)- Figure [3] shows that the 
shapes of the potential surface Eyj(ty,z), corresponding to 
the Kekule structure of the (8,8) nanotube and to the struc- 
ture of metallic phase are essentially different. For a case 
of the structure corresponding to the metallic phase of the 
(8,8) nanotube, all the barriers between adjacent minima of 
the potential surface £\y(<P,z) are equivalent. In this case, 
the same barrier AE^ = AE r should be overcome for diffu- 
sion of the fullerenes along the nanotube axis and for rota- 
tion of the fullerenes about this axis (see Figure fjp). For a 
case of the Kekule structure, two different barriers between 
the adjacent minima exist: the barrier AE^ for diffusion of 
the fullerenes along the nanotube axis, and the barrier AE r to 
rotation of the fullerenes about this axis (see Figure^). For 
the Fe@C2o fullerene, the shapes of the potential surfaces 
£w(9,z) are qualitatively the same as for C20 for both the 
considered structures of the nanotube. The calculated val- 
ues of the barriers AE& and AE r are listed in Table [T]for both 
fullerenes. The dependences of the interaction energy be- 
tween the C20 and Fe@C2o fullerenes and the nanotube on 
the relative displacement of the fullerene along the axis of 
the nanotube and the angle of its relative rotation about this 
axis, corresponding to both the considered structures of the 
(8,8) nanotube are compared as shown in Figure |4] Figure |4] 
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Table 1 - Calculated characteristics of the dynamical behaviour 
of the C20 and Fe@C2o fullerenes inside the (8,8) nanotube of 
different structure 



Nanotube structure 


Kekule structure 


Metallic phase structure 


Fullerene 


C20 


Fe@C 2 o 


C20 


Fe@C 2 o 


A£ d (meV) 


1.68 


1.73 


0.87 


0.44 


A£ r (meV) 


0.33 


0.05 


0.87 


0.44 


v d (GHz) 


72.2 


54.8 


75.7 


51.3 


V z (THz) 


0.439 


0.144 


0.558 


0.329 


V x (THz) 


1.30 


1.25 


1.36 


1.13 


V y (THz) 


2.15 


2.82 


2.08 


2.14 



AEd is the barrier for diffusion of the fullerenes along the nanotube axis, 
AE, is the barrier for rotation of the fullerenes about the nanotube axis, Vd is 
the frequency of small relative vibrations of the fullerenes along the nano- 
tube axis, V- is the frequency of small relative rotational vibrations of the 
fullerenes about the nanotube axis, v v and v v are the frequencies of small 
relative rotational vibrations of the fullerenes about the two mutually per- 
pendicular lateral axes. 

is a vivid illustration of the significant changes of the barriers 
AEa and AE r . The most dramatic change corresponds to the 
rotation of the Fe@C2o fullerene about the nanotube axis. 

Let us discuss the reason of the significant changes of the 
barriers AE^ and AE r at the Peierls transition by the example 
of the barrier AE V for rotation of the Fe@C2o endofullerene 
about the nanotube axis. Figure [5]presents the dependences 
£\v,-((p) of the interaction energies between the endofullerene 
and individual atoms of the nanotube on the angle cp of rota- 
tion of the fullerene. Figure [5] shows that maxima of depen- 
dences Ewi(tp) for individual atoms of the nanotube occur at 
different angles <p mi - and so the dependence £\v(<p) of total 
energy on the angle of rotation is essentially smoothed. In 
other words, the barrier AE r in the dependence £\y(9) of the 
total interaction energy between the endofullerene and the 
nanotube is less by an order of magnitude than the barriers 
AE T j in the dependences of the interaction energy between 
the endofullerene and only one of the nanotube atoms. Thus, 
the barrier AE V is very sensitive to the values of the barri- 
ers AE r j and angles (p m( . Therefore, the barrier AE T changes 
considerably at the found changes of the barriers AE T j and 
angles <p mi obtained from the Peierls distortions of the nano- 
tube structure and the change of the nanotube symmetry at 
the transition. It should also be noted that small barriers to 
relative motion of nanoobjects resulting from the compensa- 
tion of contributions of individual atoms to the barriers is 
a phenomenon well studied by the examples of such sys- 
tems as double-shell carbon nanoparticles lfT3l[T4l . double- 
walled carbon nanotubes Il25ll27l432l and a graphene flake 
in a graphite surface [33"]. The considerable changes of bar- 
riers for relative rotation of shells at small changes of shell 
structure were found also for double-shell carbon nanoparti- 
cles lfT3l[T4l. 

The frequencies of small vibrations of the fullerenes along 
the nanotube axis (Vd), rotational vibrations about the nano- 
tube axis (v 2 ) and rotational vibrations about two mutually 



perpendicular lateral axes (v x , v y ) are also calculated and 
listed in Table Q] The most remarkable change of frequency 
as a result of the structural phase transition corresponds to 
rotational vibrations of the Fe@C2o fullerene about the nano- 
tube axis (this agrees with the changes of the barriers). 

4 Dynamical behaviour of molecules in- 
side nanotube 

Let us consider the possible changes of the dynamical be- 
haviour of the C20 and Fe@C2o fullerenes inside the (8,8) 
nanotube caused by the structural phase transition. The 
Peierls instability transition temperature Tp was estimated 
for the (5,5) nanotube to correspond to temperature range 
T P ~ 1-15 K HEld. Both barriers AE d and AE r and the 
thermal energy kpTp are of the same order of magnitude at 
the structural Peierls phase transition (see Table [T). There- 
fore, dramatic changes of the diffusion and drift over these 
barriers can take place at the Peierls transition for the consid- 
ered pairs of the encapsulated molecules and the nanotube. 

Recently, the diffusion and drift in the periodic potential 
surface of the interaction energy dependent on the displace- 
ment z along the nanotube axis and the angle cp of the rota- 
tion about the nanotube axis were considered for the thread- 
like relative motion of walls of double-walled carbon nano- 
tubes [28 29]. The expressions for the diffusion coefficient 
and mobility of a movable wall were obtained [28 29] . In 
this study, we consider specific cases of diffusion of the 
molecules along the nanotube axis and rotational diffusion 
about this axis. In this case, the expressions for diffusion co- 
efficients, £>d and D r , mentioned above and corresponding to 
the diffusion along the nanotube axis and rotational diffusion 
about this axis, respectively, take the form: 



2 kaTj' 

Q.$ ( AEA 
A-=^exp(^— j, (2) 

where £2d and £2 r are the pre-exponential multipliers in 
the Arrhenius formula for the frequency of jumps of the 
molecule between two neighbouring global minima of the 
potential surface £\y(tp,z), 8d is the distance between neigh- 
bouring global minima for the motion of the molecule along 
the nanotube axis, 8 r is the angle between neighbouring 
global minima corresponding to the molecule rotation about 
the nanotube axis and kp is the Boltzmann constant. The 
mobility Z?d for the motion along the axis can be easily ob- 
tained from the diffusion coefficient D d using the Einstein 
ratio D&/B& = kpT. Figure [3] shows that 8d = 0.123 nm and 
8d = 0.37 nm for the (8,8) nanotube with the structure of 
the metallic phase and the Kekule structure, respectively, and 
8 r = 22.5° for the both structures of this nanotube. 

The value of the pre-exponential multiplier £2 in the Ar- 
rhenius formula is usually considered to be related with the 
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Figure 4 - Interaction energy between the C20 fullerene or Fe@C2o endofullerene and (8,8) nanotube. Dependences of the interaction 
energy £\v between the C20 fullerene (a, b) or Fe@C2o endofullerene (c, d) and the (8,8) nanotube on the fullerene displacement z along the 
nanotube axis (a, c) and on the angle (p of the fullerene rotation about the nanotube axis (b, d). Solid lines denote nanotube with the Kekule 
structure, dashed lines denote nanotube with the structure of metallic phase. The energy minimum is positioned at £\v = 0, z = and (p = 0. 
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Figure 5 - Interaction energy between the Fe@C2o endofullerene and (8,8) nanotube. Dependences of the interaction energy £w; 
between the Fe@C2o endofullerene and individual atoms of the (8,8) nanotube on the angle (p of rotation of the endofullerene about the 
nanotube axis are denoted by the thin lines. Dependence of the total interaction energy £\v between the Fe@C20 endofullerene and the (8,8) 
nanotube on the angle (p is denoted by the thick line, (a) The (8,8) carbon nanotube with the Kekule structure; (b) the (8,8) carbon nanotube 
with the structure of metallic phase. All energies are given relative to the energy minima. Only dependences £\y; with high values of the 
barriers AE r j are shown. 
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frequency V of corresponding vibrations. We suppose that 
the ratio Gl/v remains the same for relative motion of differ- 
ent carbon nanoobjects with graphene-like structure (nano- 
tube walls and fullerenes). For reorientation of the fullerenes 
of the C6o@C24o nanoparticle, the frequency multiplier Q. 
was estimated by molecular dynamics simulations having the 
value of 650 ± 350 GHz US El- We expand the potential 
surface of the intershell interaction energy near the minimum 
using the same empirical potential as in 01311141 . and calcu- 
late the frequencies of small relative librations of the shells. 
The calculated libration frequency has the value v 50 GHz, 
an order of magnitude less than that of the frequency multi- 
plier i2. In the estimations of this study, we use the values 
Cld « lOVd and Q. T w 10v z for the pre-exponential multipli- 
ers. 

The temperature dependences of the diffusion coefficients, 
D& and Z) r , estimated using expressions (f2|i are shown in Fig- 
ure [6] The dependence of the interaction energy £\v (z) on 
the displacement of the fullerene C20 has two different barri- 
ers between the neighbour minima for the case of the nano- 
tube with the Kekule structure (see Figure^). Since the fre- 
quency of the jumps of the fullerene between the neighbour 
minima exponentially depends on the barrier, the contribu- 
tion of jumps over the lower barrier into the total diffusion 
coefficient is disregarded in these estimations. The tempera- 
ture range corresponding to the Peierls transition temperature 
estimates Tp ~ 1-15 K 14,5,8] and the cases AE^/k^T < 1 
and AEr/kaT < 1, where the Arrhenius formula is adequate, 
are considered. (For rotation of the Fe@C20 inside the nano- 
tube with the Kekule structure the Arrhenius formula is not 
applicable (AE r / k%T > 1 ) at this temperature range; this case 
is considered below.) Figure [6]shows that the changes of the 
diffusion coefficients, Da and D r , at the Peierls transition can 
be of orders of magnitude. It is of interest that the diffusion 
coefficient D r for rotational diffusion of the C20 fullerene de- 
creases at the Peierls transition with the increase of tempera- 
ture. 

If a molecule is encapsulated inside a nanotube without a 
structural phase transition, the jump rotational diffusion takes 
place at low temperatures, AE^/k^T > 1, and the free rota- 
tion of the molecule occurs at high temperature, AE r /k^T < 
1. Orientational melting (a loss of the orientational order 
with an increase of temperature) has a crossover behaviour 
if the structural phase transition is absent. Firstly, orien- 
tational melting was considered for two-dimensional clus- 
ters with shell structure [34-37] and later for double-shell 
carbon nanoparticles [13, 14], double-walled carbon nanotu- 
bes [25 38 39] and carbon nanotube bundles [40|. In the case 
where a molecule is encapsulated inside a nanotube with a 
structural phase transition and the barrier AE T for rotation of 
the molecule is greater for the high-temperature phase than 
for the low-temperature phase, an inverse orientational melt- 
ing (a loss of the orientational order with a decrease of tem- 
perature) is possible. In other words, the inverse orienta- 
tional melting takes place if the Peierls transition tempera- 
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Figure 6 - Dependences of the diffusion coefficients on the re- 
ciprocal of temperature 1 /T. (a, c) The diffusion coefficient D^; 
(b) the diffusion coefficient D r . (a, b) The C20 fullerene; (c) the 
Fe@C20 endofullerene. Solid lines: The (8,8) carbon nanotube with 
the Kekule structure; dashed lines: the (8,8) carbon nanotube with 
the structure of metallic phase. 



ture lies in the range AE^ < ksTp < AE^, where AE± and 
A£ r i are the barriers for molecule rotation corresponding to 
high-temperature and low-temperature phases of the nano- 
tube, respectively. For the considered molecules inside the 
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(8,8) nanotube, these temperature ranges are estimated to 
be 3.8 < T P (K) < 10 and 0.58 < 7> (K) < 5.1 for the C 20 
and Fe@C2o fullerenes, respectively (see Table [TJ. As these 
temperature ranges are in agreement with the Peierls transi- 
tion temperature estimates Tp ~ 1-15 K JH[5l[8l, we predict 
that the inverse orientational melting is possible for the sys- 
tems considered. The inverse orientational melting should 
be more prominent for the case of the Fe@C2o fullerene with 
the greater ratio of the barriers AE r h/AE r \. 

Let us discuss the possibility of observing the changes of 
the dynamical behaviour of molecules inside armchair car- 
bon nanotubes at the Peierls transition. We believe that the 
most promising method is high-resolution transmission elec- 
tron microscopy. This method was used for visualizing dy- 
namics of processes inside nanotubes, such as reactions of 
fullerene dimerization with monitoring of time-dependent 
changes in the atomic positions |41 1 and rotation of fullerene 
chains |42|. The rotational dynamics of C6o fullerenes inside 
carbon nanotube was studied also by analysing the interme- 
diate frequency mode lattice vibrations using near-infrared 
Raman spectroscopy [43 1. The orientational melting in a sin- 
gle nanoparticle may be revealed also by IR or Raman study 
of the temperature dependence of width of spectral lines. A 
specific heat anomaly in multiwalled carbon nanotubes may 
be caused by the orientational order-disorder transition [44J. 
In the case of encapsulated magnetic molecules (for exam- 
ple, the Fe@C2o endofullerene considered above), the study 
of the temperature dependence of the electron spin resonance 
spectra could yield information on the molecule rotational 
dynamics of these molecules II4511 . 

5 Conclusive remarks 

In this Letter, we consider the changes of dynamical be- 
haviour of fullerenes encapsulated in armchair carbon nano- 
tubes caused by the Peierls transition in the nanotube struc- 
ture by the example of the C20 and Fe@C2o fullerenes inside 
the (8,8) nanotube. We apply the DFT approach to calculate 
the structure of the C20 and Fe@C2o fullerenes. The ground 
state of the (8,8) nanotube is found to be the Kekule struc- 
ture using the method of molecular orbitals. The Lennard- 
Jones potential is used for calculating the barriers for motion 
of the fullerenes along the axis and rotation about the axis of 
the (8,8) nanotube with the Kekule structure and the struc- 
ture with all equal bonds corresponding to low-temperature 
and high-temperature phases, respectively. We show that the 
changes in the coefficients of diffusion of the fullerenes along 
the nanotube axis and their rotational diffusion at the Peierls 
transition can be as much as several orders of magnitude. The 
possibility of the inverse orientational melting at the Peierls 
transition is predicted. The analogous changes of dynamical 
behaviour are also possible for other large molecules inside 
armchair nanotubes. We believe that the predicted dynamical 
phenomena can be observed using high-resolution transmis- 
sion electron microscopy, near-infrared Raman spectroscopy, 



specific heat measurements, and by study of electron spin 
resonance spectra for magnetic molecules. 

Abbreviation 

DFT, density functional theory. 
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